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Energy  Deposition  in  Microvolumes  of  Silicon 
From  High  Energy  Proton  Reactions 


1.  INTRODUCTION 

I'lic  j) I'l idiu'tiwn  nl'  I'lri  iiit  up.si'l.s  ind  lofjii.-  rrri-i'.s  by  of  liiglily  ionizing 

[i.'i  rtic'los  .such  ,<s  cosniic  ra'.  s  oc  alplia  pafticles  through  niicroelectronic  circuitry 

li.cs  boon  im.lor  uu  cstigation  fur  moi-c  tlian  five  yiuirs.  '  The  analyses  of  these 

i”.  cats  .iri'  basi'd  on  tlie  [jrobabilitv  of  a  particle  possessing  sufficient  energy  and 

stofjping  power,  and  the  probabilitv  that  the  !)article  has  a  track  segment  (chord 

lengtln  in  the  i:ollei:tion  \a)lume  long  enough  to  deposit  the  requiricl  energy  for  upset. 

The  I’hai'ges  liberated  by  the  energy  deposition  ar(.>  then  collected  as  a  pulse  or  on 

a  stcji  age  node  to  cause  a  soft  erroi-.  In  this  report,  the  ionizing  particles  are 

assumed  to  be  created  by  the  inelastic  scattei-  and  i-eactions  of  protons  in  the  energy 

2H 

range  from  100  Me\  to  1  t!e\  with  Si  nuclei.  Much  ol  the  background  for  this 
wor'k  and  some  results  are  in  the  literature  and  are  not  repeated  here. 

2.  FORMALISM  FOR  ENERGY  DEPOSITION  FROM 

NUCLEAR  REACTION  SECONDARIES 

\\  hen  describing  the  energy  deposition  by  nuclear  reaction  secondaries  in  a 
presci  ibed  convex  volume  \  ,  one  must  account  for  three  features:  (1)  the  stopping 

iKeceivcfl  for  juiblication  2a  .ianuary  |!hl4> 

I  Due  to  the  large  number  of  references  cited  ab<.)ye,  they  will  not  be  listed  here. 

•Sec  Kcfei'ences.  oaee  2a.  ) 


power-  IS  i-nefgy  dependent  ;.nd  is  variable  along  the  ti-aek;  (2>  the  i-ange  of  the 
particU-  is  finite:  id'  the  particles  are  made  locally,  that  is,  p;-oduction  occurs 
both  insute  and  outside  the  collection  volume,  and  botli  locations  contribute.  The 
fcirmalism  .issumes  tiiat  tlie  secondaries  ;ire  in  bulk  equilibrium  (the  spectrum  is 
spatially  indei>endent). 

2.1  Secondaries  Produced  Outside  the  Collection  Volume 

The  number  of  secondaries  of  a  particular  species  produced  by  incident  primary 
radiation  (for  e.\ample,  protonsi  interacting  with  the  host  nucleus  through  a  cross 
secti'in  ij  is  gi\  en  bv 

o'  =  <!>  CT  ^  P(E)  (1 ) 

2 

where  4-  is  the  primary  flux  in  protons/cm“,  a  is  the  total  cross  section  in 
(cm^/host  atom  •  proton)  for  the  particular  species,  N/A  is  the  number  of  host 
atoms/gm.  and  olbA  is  the  frequency  distribution  (probability  density)  for  second¬ 
aries  of  energy  K,  in  secondaries/Me\  .  These  dimensions  result  in  O'  being 
given  in  units  of  (secondaries /gm  •  M<-\  ). 

The  bulk  equilibrium  flux  that  results  is  given  by 

l; 

max 

/  o’(L-;  )dE 

.\'(E)  particles  iL'  ,,,\ 

- 2 -  "  - i - ; — 

cm  ■  Itle\  - j  1.  ( E)  (Me\ /cm) 

p  gm/cm 

where  p  is  the  host  density  and  I.  (E)  is  the  total  l.inear  Energy  Transfer  (LET). 

2  ' 

A  convex  surface  of  .S  cm  that  bounds  a  volume  immersed  in  an  isotropic  flux 
of  .\'(E)  particles  will  intercept 

N'(E)  S/4  particles  .  (3) 

The  probability  of  these  lying  along  a  track  (chord)  length  greater  than  or  equal 
to  the  length  required  to  deposit  ^fjj^gshold  T-ven  by  t.'(s|E,  E]).  The  function 
t'(s)  is  described  in  detail  in  Keference  6. 

Hence,  the  number  of  energy  depositions  greater  than  AE^^  from  reactions 


ic  urring  outside  the  collection  volume  is  given  by: 


max 

E  f  o(i:  »ii).  <il 

max  ^  J 

=  L,  X  - - r-rn\ - 


and  finally,  using  the  t  auchv  value  foe  mean  elioial  length,  ( 


nU  i:  )  --  /  4)  p  —  t  ts|  i;,  A  i;)) \(  I.) 

^'ai.  ^ 


where  N(E)  is  given  bv 


ME)  ^ 


/  o(E  )dE 
E 


and  has  the  normalization 


/  N'(E)dE  =  u  (mean  range) 


Thus,  the  maximum  \'alue  obtainable  by  Eq.  (5)  is 


n(Al'.  .  0)  =  4>  d  /t  4  u  . 

'  t I  :\  f 


2.2  Secondaries  Produced  inside  (he  Collection  Volume 

Tlic  numbt'r  of  seeondaries  produet’d  in  a  volume  \  by  the  same  ineidenl  flux 
of  primaries  is  given  bv 


O'  (  E)  =  <t>  Q  p  \  0(E)  pa r  tieles / MeV 


anu  the  probability  of  lying  on  an  internal  path  segment  s  of  sufficient  lengtii  to 
depositAEjj^  or  greater  is  giviui  by  (;(s|E,  AE]).  Hence,  the  number  of  energv 
di-positions  greater  than  A  E^^^  due  to  re  actions  inside  tlie  eolleetion  \olume  is  given 


n(Al;^l^)  ^  4>  a  ^  p  \  /  (dslA  E,  E))  0(|.)dE  . 


tits)  is  described  in  detail  in  Keference  t>.  Die  maximum  value  obtain aide  b\  Eq.  ' 
is  n(0)  -  4>d  4"  p  \  . 


2.3  Total  Production 


The  total  number  of  energy  depositions  greater  than  A  is  given  by  the  sum 
of  Eqs.  (5)  and  (9). 

E 


n  (AE)  =  ^  p  V  1:^  /  CisjE.  AE])N(K)dE 


^th 


E 


max 


+  /  G(s[E.AEl)  0(E)dE]  , 


(10 


The  chord  length  required  to  deposit  A  E  is  determined  from 
E-AE 


s  =  / 

E 


1 


L(E-) 


dE’i 


(11 


where  L(E-!')  is  the  electronic  stopping  power. 

In  the  lower  limit  of  the  integrals  in  n^  a  factor  g  has  been  introduced.  This 
function  is  required  for  ions  other  than  protons  because  at  low  energies 
(E  <  10  keV /nucleon)  the  electronic  stopping  power  begins  to  decline  and  atomic 
scattering  becomes  the  chief  energy  loss  mode.  Thus,  if  an  amount  of  energ> 

A  is  required  as  iotiUing  energy  deposition,  an  initial  kinetic  energy  of 
A  Ejjj/g  is  needed. 

From  the  work  of  Sattler  the  g  factor  in  silicon  is  taken  as 
g=  ^-20AE/A) 

where  A  is  the  mass  number  of  the  particle. 

2.4  DilTercnlial  Energy  Deposition  Spectra:  Event  Spectra 

The  inelastic  scatter  of  high  energy  protons  from  silicon  causes  reactions  in 
which  several  secondary  particles  are  emitted.  Hv  way  of  example  the  reaction 

Si  (p,  3p  ,  2n,  2a)  O 

involves  eight  outbound  particles  of  which  six  are  charged.  Each  of  these  has  an 
energy  deposition  spectrum.  The  particles  are  created  "simultaneously"  and 
therefore  the  total  energy  deposition  is  seen  as  a  single  event.  Thus,  to  create  tin 
correct  total  energy  deposition  spectrum  one  must  first  determini'  the  individual 


To  do  this  tine  needs 


event  spectfa  and  tlien  combine  them  in  an  appropi-iate  way. 

the  differential  of  the  n  lAK',  Kq.  (10); 

T 


-9  n_  lA  If) 
■3Al; — 


1  f 


a  ^  p  \ 


1  f(s)N(E)dL; 

-  L  il'  -  AC) 

I  Ab/g  e 


s(Ai:/g)j  xtAif/g)  ^ 


max 
/g 


(.-(slE,  AK])  0  (If)  dK 

7  L  (E  -  AE) 
e 


tl|s(AE/g){  0  (  A  E/g> 


(  13) 


wlu'i'e  we  tuive  used 
9  s  _  9  9  .s 

^  'SST. 


and 


9  s  _  1 

TaTT  ■  I.  (E  -  A  E)  • 
o 

The  separate  terms  have  been  labeled  by  Caswell'^  as  Crossers,  Stoppers, 
Starters,  and  Insiders.  C.'r'ossei-s  arc  particles  born  outside  the  collection  volume 
which  transit  coni pletel v  through  the  collection  volumes.  Crossers  are  the 
dominant  population  wherever  the  T  of  the  volume  is  small  compared  to  the  mean 
rttnge  of  the  secondary. 

Stoppers  are  particles  born  outside  the  collection  volume  whose  residual  track 
length  upon  arrival  at  the  surface  is  shorter  than  the  chord  length  along  w'hich  the 
track  lies.  Hence,  the  particle  track  terminates  in  the  interior'  of  the  collection 
volume  and  the  energy  dump  equals  the  total  kinetic  energy. 

Starters  are  particles  born  inside  the  collection  volume  whose  range  is  greater 
than  the  chord  segment  length  to  the  surface.  Hence,  they  escape  the  collection 
volume. 

Insiders  are  particles  born  inside  the  collection  volume  whose  range  is  shorter 
than  the  chord  segment  to  the  surface.  Hence,  the  particle's  entire  track  is  in  the 
interior  and  the  energy  dump  equals  the  total  kinetic  energy.  These  four  classes 
of  particles  are  shown  in  Eigure  1. 


11.  C  aswell,  11.  S.  ,  and  C'oyne,  .1.  .1.  ( 197.'))  JMicrodosimetric  spectra  and  parameters 
of  fast  neutrons,  fifth  Symposium  on  IMic rodosimetry ,  \  erbania  Hallanza, 
Italy,  J.  Hooz,  H.  (I.  Ebei't,  andH.  (1.  H.  Smith,  Eds.,  Commission  of  the 
European  C  ommunities,  Luxembourg,  El  K  54,52  d-f-c,  p.  97. 
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I'igure  1.  The  Four  Classes  of  Secondary 
Particles  Identified  by  Their  Origin  and 
Path 


}.  SUMMATION  OF  THE  INDIVIDUAL  EVENT  SPECTRA 


).  I  Convolution 

The  mathematical  procedure  for  combining  independent  spectra  is  called 
convolution.  The  independent  event  spectra  are  normalized  to  unity  and  represent 
proper  probability  densities  oi-  frequency  functions.  The  formal  rule  for  combina¬ 
tion  is 


n(6) 


(-m;, 


/ 


nj(^F)  n2(6-AE>dAE 


(14) 


where  the  n^  are  the  independent  spectra.  Caution  is  necessary  since  different 
spectra  will  usually  have  different  i-anges.  Therefore,  Eq.  (14)  must  be  divided 
into  two  parts 

6 

n(6)  =  /  n,iAE)n„(6  -AE>  dAE  (15) 

0  ' 


0  <  6  <  AE , 


(AK,  +AK„) 

1  2  max 


n(6)  =  / 

0 


nj(Ai:)  n.,(6  -  AindAK 


al:  <  6  <(ai:,  +  Ah',,) 

max  1  2  max 


Note  that  the  spectrum  with  the  least  range  takes  the  role  of  in  the  integial. 

If  there  is  a  third  particle,  a  second  convolution  must  be  performed. 

This  is  the  procedure  which  is  applied  directly  to  the  particles  of  inside  pi'oduc- 
tion.  The  starters  and  insiders  spectra  from  Kq.  (13)  are  added  and  convolved 
with  those  from  the  second  particle.  An  example  of  the  resulting  spectrum  is  shown 
in  Tigure  2. 
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f  igure  2.  linei’gv  Deposition  Spectra  af  Xeon  and  Alpha  Particle  Ueai  lion 
Products.  The  coincidencL'  spectrum  obtained  by  convolution  is  also  shown 


3.2  Weighting  of  Spectra 

The  mic  ro\  oUiine  of  the  electronic  device  that  collects  the  charge  aits  very 
much  like  a  detector.  One  can  treat  the  total  energy  event  S|)ectrum  as  the  sum  of 
se\  eral  in[>ut  sjiectra,  earn  of  these  being  identifii’d  bv  its  origin.  (  )ne  such  is  the 


-s.  V  ^  . 


•  ' '  •  - 


spectrum  that  arises  from  reactions  in  the  interior;  tiiis  spei  ti  um  is  obtaiiu  d  hv 
the  convolution  process  above.  1  rom  reactions  that  occur  outside  the  mic ruvohuiie. 
several  additional  spectra  are  created. 

In  the  example  reaction  above,  an  0*^  recoil  is  produced  vvitti  three  protons  and 
two  alpha  particles.  A  spectrum  arises  from  each  of  tfiese  plus  a  spectrum  for 
each  coincidence,  double  coincidence  and  so  fortti.  'Die  procedure  is  to  establi.'  a 
each  of  these  energy  deposition  spectra  by  the  convidution  process  and  to  assign  a 
weighting  factor  y  to  each  where  the  value  of  >  is  determined  by  the  relative  prob¬ 
ability  of  its  occurrence. 

To  determine  a  general  rule  for  assigning  a  value  of  >  one  can  proceed  as  follows: 

3.2.1  HEAVY  PA  KTICLES 

The  total  number  of  reactions  that  occur  outside  the  vidume  and  whose  heavy 
recoils  have  range  large  enough  to  reach  the  collection  volume  so  as  to  deposit 
some  energy  is  given  by 

n=$a-^Pl\  -\]  ilT) 

u  A  ‘  u  o' 

for  recoils  of  range  u,  and  where  \  is  the  volume  enclosed  b\  a  surfai  e  pa  railed 
to  V^  at  distance  u.  Tor  a  spectrum  of  particles  the  probability  of  rangi*  u  is 
r(ul  so 

u 

Y  max 

n=<t>aVP  /  -i  )  r(u)  du  .  (1 H) 

A  ^  qJ  u  o 

For  each  u  value  one  must  determine  the  probability  of  a  recoil  hitting  the 
collection  volume.  To  do  this  exactly  one  would  have  to  calculate  the  probabilities 
from  a  Monte  Carlo  geometry  (irogram.  However,  since  the  results  a  r('  to  be 
integrated  tw'ice  over,  it  is  reasonable  and  of  sufficient  accuracy  to  use  mean  values 
and  the  following  approximation. 

Consider  the  collection  volume  \  and  around  it  the  source  \olume  \  -  \  as 

o  u  o 

shown  in  Figure  3.  Between  \  and  \  create  the  parallel  surface  at  distance  x. 

“  u  o  ‘ 

The  probability  of  a  reaction  secondary  oidginating  in  this  shell  is 


.S(xl  dx 
V"  -  V 


Figure  3.  The  Source  Volume  \ 
for  any  Secondary  is  Bounded  by  a 
I’arallel  Surface  Located  at  a 
Distance  u  =  Range  Outside  the 
Collection  V  olume,  \ 

o 


PARALLEL  SURFACE  -  SOURCE  VOLUME 


From  geometric  probability  one  knows  that  the  probability  of  one  of  u  i-andom 
infinite  lines  hitting  V^.  if  the  line  hits  S(x).  is  given  by  the  ratio  of  the  surface 
areas,  S^/S(x). 

Here,  the  ranges  are  finite  and  oriented,  which  means  one  must  modify  this 
probability  by  multiplying  by  a  function  of  (x,  u)  to  account  for  the  finite  range  and 
directedness.  We  obtain 


f{x,  u)S  /2S(x) 
o 


where  the  factor  of  2  accounts  properly  for  the  directedness  of  the  track.  W  t  obtain 
by  insertion  into  Eq.  (18) 
u 


N 


max 


n  =  $  a  ^  p  /  (V,_  -  V  J  r(u)  / 


0 


u  o 


^o  ,  S(x)  dx  du 

0  u  o 


(  I'U 


The  boundary  conditions  on  f(x,  u)  are: 


(1) 

(2) 

(3) 

u 

/ 

0 


f(x,  u)  ■*  1  at  X  =  0 
f(x,  u)"*  1  as  u  *♦  <»: 

To  get  Eq.  (7)  from  Eq.  ( 

u  o 


19)  one  must  liave 

S  u/4 
o 

■  \  -  V 

u  o 


that  is. 


u 

/  f(x,  u)  dx  =  u/2  . 
0 
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Solutions  which  obey  this  condition  include 


(1)  fix,  u)  =  1  -  x/u 

(2)  =  (1  -  x/u)  +  x/u(l  -  x/u)  (1  -  2x/u) 

(3)  Power  series  in  (1  -  x/u)  x/u. 


VVe  have  used  Solutions  1  and  2  and  chosen  Solution  1  as  sufficiently  accurate. 
Solution  3  yields  nearly  identical  results  to  Solution  2.  Insertion  into  lOq.  (19) 
vields  the  nurr.ber  of  recoils  that  are  crossers  and  stoppers. 

3.  2.  2  COINCIDENCE  WITH  LIGHT  PAKTICLES 

Since  the  range  of  the  light  particles  is  so  much  greater  than  that  of  the  recoils, 
their  range  u  can  be  treated  as  an  infinite  line.  Hence  the  probability  of  light 
particle  intercept  is 

S 

o 

2snn  ■ 

The  coincidence  probabilities  of  n  particles  scoring  m  hits  are  then  binomially  dis¬ 
tributed  in  this  factor: 

{'s' 4ra’"’ '‘-ot’"""}-  <2'” 

Thus,  insertion  of  Eq,  (20)  into  Eq.  (19)  yields  the  number  of  coincidences. 

The  ratio  of  that  number  to  the  total  number  of  particles  gives  the  weighting  factor. 
Since  the  sum  of  all  the  coincidences  (including  0^*^)  equals  the  total  number  of 
particles  it  is  equivalent  to  normalize  the  weighting  by  dividing  by  their  sum.  We 
obtain 


raw 


max  u  S 


I  R(x)™  dxdu 


(21) 


The  factor  ^  cr  .  p  is  dropped  since  it  is  common  to  all  y”  .  The  factor  R(x)”’  is 
the  range  sum  distribution  and  is  the  probability  of  the  light  particle  range  being 
greater  than  x  (this  insures  the  possibility  of  coincidence]. 

The  mixing  ratios  of  weighting  factors  for  each  spectrum  are  then  determined 
by  normalizing  to  the  sum  of  the 


(22) 


y 


n 

m 


y raw 
'  m 


Thus,  in  a  reaction  which  produces  n  light  particles  and  a  heavy  recoil,  the  coin¬ 
cidence  spectrum  for  order  coincidence  (convolution)  is  weighted  by 


4.  REACTION  SPECTRA 


The  secondary  spectra  from  high  energy  proton  reactions  are  generated  by 
codes  developed  by  high  energy  physics  groups.  The  code  used  to  generate  the 
spectra  given  later  in  this  report  is  the  Medium  Energy  Collision  Cascade  Code-7 
at  Oak  Ridge  National  Laboratory.  The  work  was  performed  under  contract 
RADC-2306J320.  Other  codes  exist  which  are  suitable,  that  is,  they  produce 
reaction  secondary  energy  spectra  and  cross  sections.  The  most  notable  is  prob¬ 
ably  the  HETC  code  of  Chandler  and  Armstrong. 

The  output  of  the  codes  are  differential  spectra  in  energy  and  angle.  The 
formalism  of  chord  length  demands  isotropic  fluxes  so  we  presume  isotropic  angular 

distribution  and  the  integrated  spectra  (a  good  approximation)  over  angles. 

28 

The  output  of  a  specific  case,  a  400  MeV  proton  on  Si  ,  is  illuminating: 

(1)  There  are  70  (A,Z)  products, 

(2)  The  total  inelastic  cross  section  is  420  mB. 

(3)  The  individual  (A,  Z)  cross  sections  range  from  0.  01  mB  (several)  to 

(27,  14)  54  mB. 

(4)  There  are  evaporation  spectra  for  the  light  particle  production, 

,  .,3  „4 

p,  n,  d,  N  ,  H  . 

(5)  Thirteen  of  the  (A,  Z)  secondaries  contain  82  percent  of  the  total  cross 

section.  These  may  be  combined  by  isotopes  to  yield  eight  spectra 
sufficient  to  characterize  all  the  spectra. 

(6)  The  cross  sections  for  incident  protons  with  energies  from  100  MeV 

to  1  GeV  are  slowly  varying,  making  average  calculations 
meaningful  for  this  whole  energy  domain.  This  domain  constitutes 
the  essence  of  the  cosmic  ray  background  flux. 

The  eight  spectra  that  characterize  this  entire  secondary  population  are  listed 
in  Table  1. 
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Table  1.  Spectra  of  Secondary  {Population 


Z 

a 

a  mil 

Convolution  Product.s 

28.  27 

■1 

27 

73 

25.  26,  27 

mm 

26 

1  14 

IP 

23.  24.  25.  26 

mm 

24 

114 

2  p 

22,  23.  24 

1 1 

23 

62 

a.  p 

20.  21.  22 

10 

21 

44 

2p.  a 

16 

8 

16 

24 

2a,  2p 

14 

1 

14 

22 

2a,  3p 

12 

6 

12 

12 

30,  2p 

5.  RESULTS  AND  DISCUSSION 

Two  features  of  the  ionizine  enorfjv  deposition  in  microelertronics  are  wortiw 
of  further  mention.  Tiiese  are  (I)  the  notion  of  the  t'ollertion  volume  and  (li)  the 
charj»e  colletdion  mechani.sm  that  follows  the  deposition. 

(1)  In  the  mathematics  of  geometrical  probability  the  collection  volumes  are 
required  to  be  convex  and  are  taken  to  he  right  rectangular  volumes.  I  he  diffic.ultv 
arises  in  the  determination  of  the  size  (dimensions)  of  tlie  volume  since  several 
transport/collection  processes  are  at  w  mk.  These  are:  strong  field  collection 
within  the  depletion  zone,  diffusion  into  the  collection  volume,  and  [dasma  column 
or  "funnel"  effects.  The  chord  length-LKT  mode  of  analysis  does  not  explicitly 
treat  the  .second  two  of  the.se  proce.sse.s  and  one  can  only  adjust  the  dimensitjns  of  the 
collection  volume  in  an  a^l  hoc  W'av  to  cf>rre('t  htr  them,  lo  date  the  <letermination 
of  the  collection  volume  dimensi<.>ns  remains  an  imprecise  factor  but  thev  are 
generally  taken  as  the  lateral  flimensions  of  the  n*  colle<  tor  imfilant  (ilus  a  diffusion 
length  and  the  depletion  depth. 

(2)  The  connection  between  energy  deposition  in  the  form  of  electron -hole  pairs 

and  the  charge  required  to  represent  an  error  is  provided  througli  the 

relation: 

C^  ^  ^( coulombs)  W  (e\  c-|)  pair) 

^'■'threshold  ■  ..'  oulombs  ^ 

where  \V  -  2.  6  e\  in  silicon  and  is  the  mean  collection  efficiency  for  a  charge 
created  in  the  colUn  tion  volume.  'I  his  relationship  represents  a  noint  of  departure 
between  two  possible  modes  of  analysis: 


species.  The  diffet'ence  in  rates  can  be  traced  laryelv’  to  tlie  lit'fereine  in  cross 
section  for  speides  production.  I  he  solution  to  liq.  (10)  lor  the  proton  eva[)oration 
spectrum  is  shown  in  I'i^uta'  a.  Ni>  curve  is  shown  foi-  the  lar^jest  volume  since 
the  proton  eneryv  deposition  is  too  small  to  be  of  <  <>nsei)ucn<  e.  .As  these  i'  irves 
show,  large  numbers  of  firoton -generate. 1  etiergv  depositions  take  place  but  are  so 
small  in  value  as  not  to  .  onttdbute  to  the  soft  error  rate  (.iaslied  line). 


l  igure  .0.  IduM-gv’  Deposition 
Profiles  I'r orr.  the  Proton 
■Secondaries.  Althougli  large 
in  number,  the  firofiles  are 
concent  rate.  I  at  verv  low 
AK  values  and  so  do  not 
contribute  greatlv  to  single 
event  upset  until  256  K  HAM 
size  arravs 


The  solution  to  Eq.  (13),  the  event  spectrum  or  differential  energy  loss 
spectrum,  for  neon  is  shown  in  Figure  6.  The  two  cases  displayed,  the  largest 
and  smallest  volumes  treated,  show  the  considerable  range  of  probability  for 
energy  deposition.  In  spite  of  this  the  mean  AE  for  the  large  volume  is  onlv 
3.  6  MeV,  not  enough  to  trigger  electronic  ufjset  <-A  Ej^l^  5.  6  Ale\  ). 
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EVENT  SPECTRUM  FROM  REACTIONS 
INVOLVING  Ne7  «44mb  RECOILS 


VOLUME  3.5x14x21m 
3  6  MeV 


.  VOLUME  437SX  1.75x2.625 
\  ZE  =  .8  MsV 


0  2  4  6  6  10  12  14  16  18  20  22  24 


Figure  6.  Event  Spectra  From  Neon  Recoils  for  Two 
Volume  Sizes.  Arrows  indicate  threshold  energy 
required  for  nominal  device  upset 


As  explained  above,  the  particles  and  their  energy  deposition  can  be  identified 

by  their  origin:  outside  the  volume  or  inside.  The  relative  contribution  from  each 

source  for  the  neon  recoil  reactions  can  be  seen  in  Figure  7  for  the  largest  and 

smallest  volume.  In  the  large  volume  the  contributions  are  dominated  by  inside 

production.  This  oc:curs  because  the  mean  range  u  for  neon  is  less  than  the  mean 

c  hord  length  T:  u.,  =  .'1.6  u,  r,  =5  u.  In  the  limiting  form  the  contributions 

Ne  larger  _  ” 

of  outside  to  inside  productions  approach  u/T  to  1.  This  effect  shows  clearly  in 
the  small  box  where  outside  product’  -n  dominates.  Here  =  1.25  u. 
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WEIGHTED  CONTRIBUTIONS  TO  C  ‘  SPECTRUM 
■  ■■  RECOIL  C'^  t53-  34 

RECOIL  C'^ -I-  a  COINCIDENCE  83=41 
It «  »  RECOIL  C'*  +  20.  COINCIDENCE  bj  “  2 1 
000  RECOIL  C'*  +  3a  COINCIDENCE  b,  *  04 
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Figure  8.  Weighted  C  and  Coincidcr;  e  Kvent  Spectra,  \\  ith  Weighting 
Fractions  Determined  bv  the  Procedure  Detailed  in  the  'I'ext 


Hinallv.  the  complete  result  ;)!'  this  investigation  is  shown  in  Figure  9.  ’I  his 

is  the  grand  summation  of  all  tontributions.  Si  inelastic  recoils,  Al,  anti  Mg 

recoils,  plus  the  evaporation  particles  <v,  p  are  added  to  the  coincidence -convoked 

results  for  Na,  Ne,  O,  N,  and  C.  An  event  rate  stealing  curve  is  addeil  to  shtiw  the 

expeetetl  total  upset  rate.  The  ordinate  has  units  of  upsets  per  proton  per  cm^ 

-  12 

per  second  per  cell  of  volume  indicated.  The  intercept  value  of  2.  6  X  10  for 

5 

the  V  (4K  DHAM)  case  yields  ~9  events  per  vear  for  a  10  DRAM.  This  repre- 
o 

sents  about  5  percent  of  the  observed  event  rate,  the  remainder  di'e  to  heavy  cosmic- 

ray  ionizing  tracks.  The  plot  shows  that  this  event  rate  dec  lines  by  a  factor  of 

about  3  for  the  smallest  volume.  The  importance  of  the  scaling  law-  is  equally 

c  lear.  If  the  decrease  in  the  threshoki  energy  for  upset  does  ncit  follow  a  1/s 

(s  is  the  dimension  reduc  tion  fac  tor)  rule  but  a  I /s'^  rule  then  a  dramatic  increase 

in  the  event  rate  would  occur.  Lastly,  the  resolution  of  these  calculations  is 

0.  1  Mc\  .  Thus,  the  ultimate  domain  of  \  l-SI,  that  with  AK  ,  below  0.  1  MeV,  is 

th 

not  addressed. 


PROTON  CONTRIBUTION 
NOT  INCLUDED 


VQ=3  5jM4x2iAi 


EVENT  RATE 
SCALING  LINE 
IF  AE  -  t/S^ 


Vt=  1  75*7*11 5 


Vj=  075*3  5*5  75 


437Siit  7Si2.S7S 


i''igure  9.  Ciratul  Summation  for  Invents  Greater  Than  A  !■;  for  Hii;h  l^ner^v• 
Proton  (100  Me\'  to  1  Cie\')  Heai'tion  I'rodurts  in  Silicon.  Scalintj  line 
intercepts  v'ieid  event  rates 


6.  VALIDITY  -  EXTENSION  TO  VERY  SMALL  SIZES 
AND  VERY  LOW  THRESHOLD  ENERGIES 


From  the  outset,  this  calculation  has  been  based  on  the  notion  that  it  would  Ite 
sufficient  to  characterize  the  eneryv  deposition  and  thus  upset  rates  hv  treating  onlv 
the  major  contributions  to  the  variance  of  the  eneri>v  defjosition  dist ri hutions.  ’1  he 
cosniic  rav  work  was  dominated,  therefore,  bv  the  t-hord  length  —  1  dlT  contribution. 
In  this  work,  with  lower  enerjjy  parti<  les,  it  was  necessarv  to  iiu  lude  effects  of 
finite  range  and  changing  I.F'l'  (linear  en(>rgv  transfer)  along  the  path  in  addition. 

.As  one  proceeds  to  even  smaller  dimensions  it  will  become  necessarv  to  include 
two  more  contributions:  (1)  (ielta  rav  pro<iuctioti  with  a  track  no  longer  rc'present- 
able  bv  a  line  but  bv  a  3-1)  fuzzv  i-vlinder  with  large  internal  gradients;  and 
(2)  energy  straggling  effects. 

The  entire  domain  of  contribution  can  l>e  shown  on  a  single  diagr-am.  I  ollow  - 

1  g 

ing  the  lead  of  Kellerer  and  C'hmrdi'vskv  “  one  cati  i-reate  a  figure  like  I  igure  ID. 


12.  Kellerer.  A.  M.  .  and  C'hmelevskv.  I).  (  IDV  .o)  C'ril('i  ia  for  the  appli  cahi  I  i  t  \ 
Idi  r,  Had.  Hes.  .  6:);226. 


These  diagrams  have  neat  boundaries  only  for  spherical  shapes  so  I  retain  those. 

The  ordinate  is  size  (diameterl  and  the  abscissa  is  energy/nucleon.  This  diagram 
1 6 

is  for  0  ions.  The  domain  boundary  lines  are  attempts  to  show  where  the  eliect 
contributes  10  percent  to  the  variance  of  the  energy  deposition. 


OS  i  2  .«  1  2  6  10  20  60  100  200  600  1000 

MeV/NUCLEON 


Figure  10.  Oomain  Diagram  for  Contributions  to  tlie  \  ariance  of  the  Fvent 
Spectra.  Each  "boundarv"'  represents  '  10  percent  contribution  conipare  I 
to  linear  energv  transfer  chord  length  contribution 

The  single  event  upsets  began  in  a  doniain  of  relativelv  large  size  and  cosmic 
ray  energy,  henc''  ipper  right  on  diagram.  Scaling  and  size  reduction  proceed 
downward  verticallv.  Inclusion  of  nuclear  reac  tion  secondaries  of  low  energv  and 
finite  range  spanned  domain  11  and  into  I.  Further  size  reduction  will  l  arrv  into 
domain  III  and  I\  ,  delta  rav  and  finite  track  cross  sei-tion  distribution  and  energv 
straggling.  This  diagram  shows  that  the  formalism  of  this  work  does  not  need 
modification  for  use  in  advanced  \  LSI  size  domains  (<  Ip)  but  that  the  cosniic’ 
interactions  do  require  a  modification  to  account  for  the  .'1 -D  cliai’acleristics  of 
charge  distributions  around  the  axis  of  the  track  (delta  i-avs)  and  subsequentlv  to 
reflect  energv  straggling. 


7.  COMPUTER  CODE  EDPMON 


T  he  in  liviTuiil  eompuler  pt'OijrHms  \isi-.l  in  tliese  enlenlatiinis  ;uiT  in  'ofeeninL' 

1  aU'ulations  if  Ixcfei'eneos  I,  6,  aiici  8  havi-  been  i-ombine.l  iiuo  one  pioai’a’"  ealle.l 
I'.'nPMOX.  T  his  tr.oi  e  etTieienl  i-ode  was  written  i'v  Air.  Alark  Aloses  of  1T('  iloi'T 
Ht'seareh  unTer  LSAK  fontraet  I-' ll'628 -82 -C -0  IT!'.  IIDI'AIO.X  will  corr.  pate  on 
input  of  the  constants  of  the  problem  (1)  the  chord  lontjth  distribution  .'unctions  an  l 
first  and  second  moments  of  the  di.stributions  for  rectanyular  \'olumps.  ((trovision 
is  made  to  read  in  chord  length  table.s  ba.sed  on  other  geometries),  (2)  the  event 
rate  of  energv  defjosition  from  a  single  particle  species  is  proiluced  in  hulk 
'■tlvtiUbrivim  (so  called  integral  mode).  (2)  the  differential  energv  defjosition  spectrum 
!or  recoil  and  light  secondarv  particles,  (these  calculations  are  (performed 
scparatelv  for  particles  born  outside  the  collection  volume  and  (larticles  born 
ir.side  the  collection  volume).  (4)  coincidence  spectra  up  to  triple  coincidence  (via 
convolut  ion  I .  and  (5)  first  and  secoiid  moments  of  the  final  energv  deposition 

J  g 

spi'ctra.  -An  example  of  a  single  .speeies  run,  O  -2  0.  using  single  and  double 
.  oincidenc('.  used  206  seeonds  of  Central  Processor  time. 
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